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We study fully synchronized states in scale-free networks of chaotic logistic maps as a function of both 
dynamical and topological parameters. Three different network topologies are considered: (i) random scale-free 
topology, (ii) deterministic pseudo-fractal scale-free network, and (iii) Apollonian network. For the random 
scale-free topology we find a coupling strength threshold beyond which full synchronization is attained. This 
threshold scales as fc"'', where k is the outgoing connectivity and /i depends on the local nonlinearity. For 
deterministic scale-free networks coherence is observed only when the coupling strength is proportional to the 
neighbor connectivity. We show that the transition to coherence is of first-order and study the role of the most 
connected nodes in the collective dynamics of oscillators in scale-free networks. 
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I. INTRODUCTION AND MODEL 

Recently, intensive research on the structure and dynamics 
of networks has provided insight for many systems where they 
arise naturally QlHI^l- Complex networks appear in a wide 
variety of fields ranging from lasers 0], granular media @, 
0], quantum transport 0], colloidal suspensions 0], electrical 
circuits 0], and time series analysis llo||. to heart rhythms 

e, epidemics Il 4| l3 ] .protein folding I14ll . and locomotion 
among others ulSllll- 
From the mathematical point of view, a network is a graph, 
composed by nodes or vertices and by their connections or 
edges 0]. When studying network dynamics one frequently 
assumes a regular structure where each node evolves ac- 
cording to some more or less com plic ated dynamics, typi- 
cally fixed points fl6ll . limit cycles llTll or chaotic attractors 
El El. When studying network structure, one usually ne- 
glects node dynamics and all complexity is introduced by the 
way nodes are connected to each other, i.e. by the network 
topology. With respect to their topology, networks are usually 
divided into three large classes |j3]: random networks, where 
all the nodes are randomly connected i20ll . small-world net- 
works introduced recently by Watts and Strogatz l2ll l23l as 
a middle ground between regular and random networks, and 
the scale-free networks (e.g. Barabasi and Albert |(23j|), where 
growth and preferential attachment are considered. 

The next logical step toward real network dynamics is to 
consider simultaneously structural and dynamic complexity. 
One important question addressed in this context is to know if 
synchronization between oscillators in such complex topolo- 
gies would appear and under which conditions it prevails. In 
fact, coherent behavior of oscillator networks with complex 



topologies has been studied for the random topology 124125 
and small-world topolo gv 12 ^ ii27| |. i28i. iSl . However, apart 
from a few exceptions 1251 l30l 13111 . there is a quite general 
lack of studies tackling synchronization of chaotic oscillators 
in scale-free topologies. 

In this paper we present detailed results concerning syn- 
chronization in oscillator networks with scale-free topologies. 
Our purpose is to determine under which conditions scale-free 
topologies enable the emergence of coherent behavior As a 
general result, we present evidences that the transition to syn- 
chronization is of first-order. Our model reads 



(1) 



where i = 1, . . . , L and t label discrete space and time respec- 
tively, L being the total number of oscillators, < e < 1 is 
the coupling parameter, ICi represents the set of labels of the 
neighbors of node i, ki represents the number of such neigh- 
bors, and Mi = SjeK; ^i* normalizes the interaction term in 
Eq. Q. The function f{x) is a continuous function governing 
node dynamics when connections are absent. Here we choose 
the well-known quadratic map f{x) = 1 — ax^, where the 
free parameter a is restricted to the interval —0.25 < a < 2 
and contains all possible dynamical regimes from a fixed point 
(e.g. a=0) to fully developed chaotic orbits (e.g. a = 2). The 
parameter a is a real number controlling the homogeneity 
in the coupling: positive values of a enhance the coupling 
strength with sites having larger number of neighbors, while 
negative values favor sites having less neighbors. For a = 
the coupling between each site and its neighborhood is homo- 
geneous, i.e. it is independent on the coordination. 

The linear stability of the coherent states xt,i = X, Vi is 
governed by the variational equations of Eq. Q, whose diag- 
onal form reads Isllsllsl 
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where A{e\i) is the Lyapunov exponent, Df{X) represents 
the identity matrix muhipHed by the derivative of f{x) com- 
puted atx — X and are eigenvalues of the coupling matrix 
G whose diagonal values are gu = 1, while off-diagonal el- 
ements are gij — —k'^'/Mi if nodes i and j are coupled and 
zero otherwise. If G has zero-sum rows and all its eigenval- 
ues Ai < A2 < • ■ • < Al are real, then Ai = corresponds 
to the mode parallel to the synchronization manifold and the 
largest Lyapunov exponent defines a master stability function 
II32I1 . The coherent state is stable whenever A(eAi) < for 
i = 2, . . . , L ll32l I33I I34I1 . In our case, it is easy to check 
that indeed G have zero-row sum, yielding Ai = and all its 
eigenvalues are real, since det(G — AI) = det((G — AI) where 
G is a symmetric matrix, namely G = HKLKH with L being 
the Laplacian of the network ll32ll35ll . and matrices H and K 
being the diagonal matrices with elements Ha = TV; and 
Kii = fc"^^ respectively. 

From Eq. (|5J, taking into account the ordering of the eigen- 
values A;, one easily concludes that the stability condition for 
chaotic maps reads 

^-7('^^ <.< ^ + 7^'^\ (3) 
A2 Ai 

where A is the Lyapunov exponent of the local map. In partic- 
ular there is a range of coupling strengths enabling synchro- 
nizability if < (1 +e"^)/(l -e"^) holds. Therefore, 

by computing the eigenvalues of the Laplacian matrix one is 
able to find the range of couplings for which coherent states 
are stable. For more detailed results see Ref. 1361 . 

Instead of starting from coherent states and studying their 
stability to perturbations, in this paper we consider large sam- 
ples of random initial configurations and study how much and 
under what conditions they converge toward a coherent state. 
This procedure not only reveals the existence of stable solu- 
tions but also gives a rough measure of its basin of attraction. 

Earlier results ll25ll show a transition to full synchronization 
for two particular values of the nonlinearity a in the homo- 
geneous regime (a = 0), when either the coupling strength 
or the number of outgoing connections are varied. Here, we 
show that the threshold value of such a transition as a func- 
tion of coupling strength and outgoing connectivity obeys a 
power-law with an exponent that depends on the nonlinear- 
ity. We study not only the usual random scale-free network 
of Barabasi and Albert ll23ll . but also deterministic scale-free 
networks constructed in an iterative way i37l l38l I39I1 . De- 
terministic scale-free networks are analytically easier to han- 
dle iIstI I40I1 . Deterministic networks are applied for instance 
in^^n^ systems l39il . and geographical and social networks 

We consider in Section |lll the homogeneous coupling 
regime for the random scale-free topology. In SectionHIIlwe 
extend our results to two deterministic scale-free networks, 
namely a pseudo-fractal network l38ll and an Apollonian net- 
work [39|. Discussion and conclusions are given in Section 




FIG. 1: Typical histogram of the standard mean square amplitude 
deviation as a function of the coupling strength e, showing the 
transition to coherence for a sample of 500 initial random config- 
urations. Here A'^ represents the fraction of configurations, and we 
discarded transients of 10* time steps and fixed nonlinearity a = 2, 
connectivity k — mo = 8, and number of nodes L = 1000, and 
a = 0. The base of the logarithm is 10. 
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FIG. 2: Typical histogram of the standard mean square-deviation 
logj^Q as a function of connectivity k, for e — 0.95, a = 2, a = 
and L — 1000. The same conditions and initial configurations of 
Fig.Qwere used. 



II. RANDOM SCALE-FREE NETWORKS 

Random scale-free networks share with many real net- 
works, e.g. the World Wide Web, two generic mechanisms: 
growth and preferential attachment 0]. In this Section, we 
use the algorithm of Barabasi and Albert 0, Esll to construct 
the network: starting with a small number of nodes, say m.Q, 
fully interconnected with each other, one adds iteratively a 
new node with k new edges, which connect randomly the new 
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node with previous nodes, depending on their own number of 
connections. As a general feature, one finds j^] a connectiv- 
ity distribution which follows a power law with an exponent 
7 = 3, independently on mo and k. After a certain number of 
iterations, ones has a network with L nodes, and then we place 
chaotic maps at the nodes, according to Eq. Q, and observe 
if they synchronize or not after some transient. 

A suitable approach to study synchronization of chaotic os- 
cillators on an arbitrary network topology lEsIl is to compute 
the standard mean square-deviation 
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where Xt is the average amplitude at a given time step t. As 
one easily sees, all the nodes are synchronized at the same 
amplitude whenever is zero within numerical precision, 
i.e. (T^ ^ 10~'^°. We call these fully synchronized states co- 
herent states, to distinguish them from partially synchronized 
configurations, when several different clusters of nodes with 
the same amplitude are observed O- 

In a previous work by Jost and Joy l25ll . concerning lattices 
of coupled maps with different coupling topologies, a tran- 
sition to coherence between chaotic maps was found, when 
considering the Barabasi-Albert network, occurring for par- 
ticularly high coupling strengths, typically of the order of 
£c ^ 0.9. Our simulations have shown that these transitions 
occur after discarding transients of ~ 10''^ time steps and they 
do not change significantly with the network size. Moreover, 
this transition to coherence is robust with respect to initial con- 
figurations. 

Figure [2 shows a typical histogram of the standard mean 
square-deviation as a function of the coupling strength e, com- 
puted from a sample of 500 initial configurations, and fixing 
L = 1000, a = 2, and k = mg = 8. From the histogram, 
one clearly sees the sharp transition to coherence and also its 
robustness to initial configurations, since for each coupling 
strength all the final configurations have approximately the 
same standard mean square-deviation. In particular, above 
the threshold Sc ^ 0.9, all initial configurations converge to- 
ward a coherent state, indicating that in this parameter region 
the basin of attraction of coherent states fills almost the entire 
phase space. These two features, sharp transition to coher- 
ence and robustness with respect to initial configurations, are 
also observed when varying the connectivity fc, as illustrated 
by Fig.|2] Both Figs. ^ and |2] were drawn fixing one of the 
parameters, e or k. Our simulations show that for the fully 
chaotic regime (a — 2) the transition to coherence occurs for 
gradually smaller coupling strength if the connectivity k is in- 
creased. Figure |3t illustrates this fact, plotting the fraction 
Na-=o of initial configurations which converge to a coherent 
state. One sees a clear transition to coherence. Computing 
similar histograms for other values of a, smaller then a = 2, 
and projecting them in the (e, k) plane one observes similar 
transition lines in ranges with smaller coupling strengths. Fig- 
ure |3j) illustrates this fact by plotting the threshold values, 
Ec and kc, at the transition curves for (from bottom to top) 
a = 1.5, 1.6, 1.8, 1.7, 1.9 and 2, in the same conditions as in 
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FIG. 3: Transition to coherence as a function of connectivity k 
and coupling strength e. (a) Fraction Ncr^o of coherent states 
from 500 random initial configurations for a = 2. (b) Coher- 
ence transition curves in the (e, k) plane for (from bottom to top) 
a = 1.5, 1.6, 1.8, 1.7, 1.9 and a — 2, and (c) the same transition 
in a log-log plot (base 10), showing power-law dependence between 
connectivity and coupling strength for the transition curves, with an 
exponent /j, which depends on the value of a (see inset). Here a = 0, 
L = 1000 and we used transients of 10* time steps. By increasing 
the transient size to ~ 10® one sees clearly that the transition to co- 
herence is of first-order either (d) when varying the coupling strength 
£ or (e) when varying the outgoing connectivity k. 



Fig. |3p. For all these values of a, the single uncoupled map 
shows chaotic orbits, or at least the orbits have very large pe- 
riods T > 10^. Note that the curve for a = 1.8 is below 
that for a = 1.7; this slight discrepancy is due to the fact that 
for a ~ 1.8 the amplitudes of the logistic map vary (chaot- 
ically) in a smaller interval than that observed for a = 1.7. 
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As illustrated in Fig.|3j;, all curves obey, within our statistical 
precision, a power-law, 

Sc oc k-^\ (5) 

For the six values of a above, the exponents are respectively 
/i = 0.2345,0.2354,0.2353,0.2231,0.2023 and 0.1804. In 
other words, the exponent is almost constant below a ~ 1.7 
and decreases above this value, as illustrated in the inset of 
Fig.ia. 

In order to determine the nature of the transition to coher- 
ence seen in Fig.|3^, we show in Fig.lSjl a high-resolution plot 
of A^cr=o as a function of e for different connectivities. Here 
one clearly sees a well-defined jump indicating that the tran- 
sition to coherence is of first-order One also observes first- 
order phase transitions when the outgoing connectivity k is 
varied (see Fig.|3^). That transitions are indeed of first order 
is easily recognized by the clear existence of hysteresis: when 
increasing either e or fc the configuration eventually falls into 
a coherent state, no longer spontaneously desynchronizing, no 
matter how far the parameters are tuned back. 

In this Section we only consider the case of homogeneous 
coupling (a = 0). For a > 0, when the coupling to nodes with 
large number of neighbors is strengthened, we find transitions 
to coherence similar to the ones illustrated in Fig.|3] occurring 
at weaker coupling strengths. 

As a general conclusion one could say that, although the ex- 
ponent 7 of the power-law distribution of connections charac- 
terizing scale-free networks does not depend on the outgoing 
connectivity k synchronization behavior is quite sensi- 
tive to this quantity. 

III. DETERMINISTIC SCALE-FREE NETWORKS 

In the previous Section we focused on random scale-free 
networks, i.e. growing networks where new nodes are con- 
nected following probabilistic rules. Although this stochastic- 
ity is typical for real networks, it is more difficult to handle 
analytically IjtIi . In this Section we study a different type of 
networks; deterministic scale-free networks |37, 38, 39j. 

In particular, we use two different deterministic topologies, 
namely the pseudo-fractal scale-free network introduced by 
Dorogovtsev et al Issll . which is similar to the first determin- 
istic scale-free network proposed by Barabasi et al IstIi . and 
was recently applied, e.g. to studies of opinion formation l4l|l. 
and the Apollonian network introduced by Andrade et al 13^ . 

The pseudo-fractal network of Dorogovtsev is obtained, 
starting from three nodes interconnected with each other, and 
at each iteration each edge generates a new node, attached to 
its two vertices. Figure]^ illustrates this network after three 
iterations, i.e. with three generations of nodes. With such a 
construction the number of nodes L„ and the number of con- 
nections Mn increases as l38ll 

Ln = 1(3" + 1) , (6a) 

M„ = 3"+i , (6b) 

where n is the number of iteration steps (generations). More- 
over, at iteration n the number of nodes with degree k = 




FIG. 4: Illustrations of two deterministic scale-free networks: (a) the 
pseudo-fractal network and (b) the Apollonian network 13^ . 
Identical symbols label nodes belonging to the same generation n 
(see text), namely Q for ?i = 0, ■ for n — 1 and • for ?i = 2. 



2, 2^, . . . , 2"-\ 2" and 2"+^ is equal to 3", 3""\ . . . , 3^, 3 
and 3 respectively, yielding a power-law distribution with ex- 
ponent 7 = 1 + In3/ln2 c^i 2.585. 

The Apollonian network has a construction algorithm dif- 
ferent from the pseudo-fractal network: one starts with three 
interconnected nodes, defining a triangle. At n = one puts 
a new node at the center of the triangle, joined to the three 
other nodes, and thus defining three new smaller triangles. At 
iteration n ~ 1 one adds at the center of each of these three 
triangles a new node, connected to the three vertices of the 
triangle, defining nine new triangles; at iteration n ~ 2 one 
adds one new node at the center of each of these nine trian- 
gles, and so on (see Fig.|4j5). With this construction procedure 
one obtains a deterministic scale-free network C9II . where the 
number of nodes L„ and the number of connections M„ are 
given respectively by 

Ln = i(3"+i + 5), (7a) 
M„ = |(3"+i + l). (7b) 

At iteration n, the number of nodes with degree k = 
3,3 • 2,3 • 22,...,3 • 2"-i,3 • 2" and 2"+^ is equal to 
3", 3"~^, 3"~^, . . . , 3^, 3, 1 and 3 respectively, yielding a 
power-law distribution with the same exponent 7 as the one 
found for the pseudo-fractal network. 

From Fig.l^and the description above, one easily concludes 
that for the pseudo-fractal network the outgoing connectiv- 
ity is fixed at fc = 2, while for Apollonian networks one has 
k = i. Despite both networks have a small number of outgo- 
ing connections, they are quite different from the geometrical 
point of view. In fact, while the pseudo-fractal network has 
no metric, Apollonian networks are embedded in Euclidean 
space and fill it densely as n cx), b eing particularly suited 
for describing geographical situations 13911 . 

As mentioned by Barabasi et al a strong advantage 
of deterministic networks is that it is often possible to com- 
pute analytically their properties, for example the adjacency 
matrix, whose eigenvalue spectrum characterizes the topol- 
ogy 1121. A simple way to write the adjacency matrix of the 
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pseudo-fractal network is 



Ai-1 Mn-l 



(8) 



L„xL„ 



where Ln is given by Eq. ( l6al . represents the transpose 
matrix of Ai and for each generation n = 1,2,... the matrix 
7VJ„ reads 



Mr 



Mn-l Mn-l 



Bn-1 



(9) 



with 
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A ... 
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and whose starting form is 



1 1 

1 1 
1 1 



(11) 



3x3 



For the Apollonian network, the adjacency matrix is given 
by the same recurrence of Eq. but this time with 



111 
10 11 
110 1 
1110 



(12) 



and Mn being a matrix with (3" + 5)/2 rows and 3" columns 
and having in each column three non-zero elements only. 

For stability analysis purposes (see Section|IJ, one could de- 
rive the Laplacian matrices directly from these adjacency ma- 
trices, multiplying the adjacency matrices by —1 and adding 
the appropriate number of connections of each node i along 
the main diagonal. 

As shown in Fig. |5] despite having quite similar struc- 
tural properties 11381 13911 . the global dynamics of the entirely 
deterministic scale-free networks shows quite different be- 
havior than the one observed for the Barabasi-Albert model 
in the previous section. Namely, there is no coherence ob- 
served for the fully chaotic map for a = 2. In fact, from 
Fig-Elone sees that the standard mean square-deviation never 
vanishes. Instead, it is characterized by some large value 
which is almost constant beyond the weak coupling regime 
(e > 0.2). In the weak coupling regime (e < 0.2) the standard 
mean square-deviation is even larger, since the coupling is not 
strong enough to compensate the highly chaotic local dynam- 
ics (a = 2). Our simulations have shown that this feature re- 
mains valid for any transient up to 10^ time steps, and it seems 
to be valid for any value of a for which the quadratic map 
supports chaotic orbits. One possible physical explanation for 
this absence of synchronizability is that long range random 
connections are crucial to improve the ability for synchroniza- 
tion and, due to the deterministic construction of the network. 



there are no long range connections as in the Barabasi-Albert 
scale-free network. Rigorously speaking, the absence of syn- 
chronizability is valid only within the range < e < 1. Of 
course that, as long as condition (|3} holds, there is for sure 
some range of coupling strengths for which the coherent states 
are stable. However, since we are working with maps of the 
interval, we neglect coupling strengths outside the unit inter- 
val, otherwise it is not possible to guarantee convergence for 
all initial configurations. 

We plot the average (logj^Q cr^) as a function of both the 
nonlinearity a and the coupling strength e simultaneously. 
Figure|S]shows this dependence for the full range of the cou- 
pling strength < £ < 1 and for nonlinearities above the 
accumulation point of the first period-doubling cascade of the 
quadratic map, namely a ~ 1.411. InFig.|6^ we use a pseudo- 
fractal network, while in Fig. |5J) the Apollonian network is 
considered. In both cases six generations of nodes are taken, 
yielding a total of L = 1095 nodes for the pseudo-fractal net- 
work and L = 1096 for the Apollonian network. 

As can be seen from these figures, in both cases one has two 
main regions: (I) a region where the standard mean square- 
deviation is large and varies smoothly with the parameters and 
(II) a region where the mean square-deviation is smaller but 
has larger fluctuations. For the Apollonian network the irreg- 
ular region is characterized by significantly smaller values for 
the standard mean square-deviation. 

The results observed in the histograms of Fig.|6lare some- 
how surprising, since irregular variations of the standard mean 
square-deviation occur for low nonlinearity and high coupling 




FIG. 5: Typical histogram of the standard mean square-deviation 
logjQ cr^ for the pseudo-fractal network as a function of the cou- 
pling strength e. Similar result is obtained for the Apollonian net- 
work. Nodes are ruled by the map f{x) = 1 — 2x^ and we fixed 
a = 0. For the pseudo-fractal network we fixed the number of 
nodes L = 1095 (6 generations of nodes), while for the Apollo- 
nian network L = 1096 (6 generations of nodes). For each coupling 
strength, the value A'^ indicates the fraction of 500 initial configura- 
tions, and we used transients of 10^ time steps. 
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FIG. 6: Histogram of the standard mean square-deviation as a function of nonlinearity a and coupling strength! e, for deterministic scale-free 
networks, namely (a) pseudo-fractal network and (b) Apollonian network. The mean square-deviation is averaged over a sample of 500 initial 
configurations and during 100 time steps, after discarding transients of 10^ time steps. Here q = and the base of the logarithm is 10. 
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FIG. 7: Transitions to coherence in deterministic scale-free networks, when synchronizing the first g generations of nodes out of I generations 
(see text), (a) pseudo-fractal network and (b) Apollonian network. The collective dynamical behavior is quite insensitive to hubs (see text). 
Insets show that transitions to coherence are of first-order. For each network, we use = 9 generations of nodes and a = 2 fixed. The base of 
the logarithm is 10. 



strengths, precisely where one would expect the most regular 
behavior of the node dynamics. However, this irregularity is 
just apparent, since (logj^Q a^) is an average over a sample of 
initial conditions. Whenever some initial configuration leads 
to coherence the zero standard deviation decreases this aver- 
age. Therefore, for a < 1.7, i.e. in the region of irregular vari- 
ations of (T^, coherent solutions are observed. In fact, from the 
stability condition in Eq. Q one sees that for periodic maps, 
the lower boundary of the e-range is always negative while 



the upper is positive, yielding always a finite range of cou- 
pling strengths where synchronizability is possible. Since for 
Apollonian networks, the fluctuations occur at small values 
of the standard deviation, this means that there is probably a 
larger number of coherent solutions. 

Figure |6l indicates that there is a lack of coherent solutions 
above a ~ 1.7. To explain this fact one should notice that both 
the pseudo-fractal and Apollonian networks have small outgo- 
ing connectivities, k — 2 and fc = 3 respectively. Since one 
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FIG. 8: Inducing transition to coherence by varying the heterogeneity 
a (see Eq. Q) in scale-free networks, (a) pseudo-fractal network and 
(b) Apollonian network. For strong heterogeneity coherence appears 
beyond a relatively high coupling strength, but it disappears again 
for very large couplings (see text). For each network, we use £ — 6 
generations of nodes and fix a = 2. (c) and (d) show high-resolution 
plots of as a function of e for a = 2, emphasizing the first-order 
phase transition to coherence. The base of the logarithm is 10. 



also observes almost no coherent solution for random scale- 
free networks neither for k ^ 2 nor for k ~ 3 (see Section 
IIH . we believe that the outgoing connectivity k is the main 
parameter controlling synchronization between oscillators in 
complex networks. By choosing another deterministic scale- 
free network with a higher outgoing connectivity, say k = 10, 
one might see coherent solutions beyond a coupling thresh- 
old value approximately similar to those computed for random 
scale-free networks (see Fig.|3jl. 

In the remainder of this section we will consider the two 
deterministic scale-free topologies, and study possible ways 
of inducing coherent states. 

Starting from a total number of £ generations, one efficient 
way of inducing coherence is by imposing synchronization 
among a certain number of g < £ generations. By generation 
we mean the set of new nodes appearing simultaneously at 
a given iteration n, during the 'construction' of the network. 



For instance, in the Apollonian network, the first generation 
has Li = 3 nodes, the second has L2 — Li = 9 and the nth 
generation has L„ — L„_i = 3" nodes. In other words, we are 
interested in the collective effects when the most connected 
nodes (hubs) are synchronized. We will show that, in fact, 
hubs play no dominant role for the synchronization of nodes. 

Figure shows the standard mean square-deviation as a 
function of coupling strength for pseudo-fractal (Fig.^) and 
Apollonian networks (Fig. 03). In each case we choose the 
fully chaotic map (a = 2) and impose synchronization among 
the nodes of the first g generations by setting them to be their 
mean amplitude at each time-step. 

For both networks, one sees from Fig. that the standard 
mean square-deviation remains large when synchronization is 
imposed to all g < ^ — 2 generations. Coherent solutions are 
only observed for g = t — 2 and g = t~l, beyond a coupling 
threshold which is smaller for the latter case. Surprisingly, 
for g = t — \ the transition to coherence occurs precisely for 
the same coupling strength in both networks. This may be 
due to the fact that, the fraction Lg/L( of nodes on which one 
imposes synchronization is approximately the same for both 
networks. For g = £ — 2 the pseudo-fractal network shows co- 
herence only above very high coupling strengths, near e ^ 1, 
while for Apollonian networks the threshold is much lower. 
Although in this case the fraction of nodes on which one im- 
poses synchronization is also similar for both networks, it is 
much smaller than in the case where g = £ — 1. The transition 
to coherence occurs at different coupling strengths because 
the number of synchronized nodes is not enough to suppress 
the effect of the outgoing connectivity. So, since the outgoing 
connectivity is larger for the Apollonian network, its transition 
to coherence occurs for weaker coupling strengths. For both 
networks, one obtains similar results for any higher value £ of 
generations since the quotient of the number of nodes between 
two successive generations L„/L„_i — > 3 as n increases. 

As a general remark, one observes from Fig. that one 
needs to synchronize a rather high number of generations to 
induce coherence. Therefore, it seems that, dynamical col- 
lective behavior on scale-free networks is quite insensitive to 
hubs. As shown in the insets of Fig. ^ and03, the transition 
to coherence is of first-order. 

Another way to induce coherence in these two deterministic 
scale-free networks is, instead of imposing synchronization to 
the most connected nodes, to strengthen their coupling to the 
other nodes by taking a > in Eq. Q. Figure|8lillustrates the 
transition to coherence by varying the heterogeneity a for the 
pseudo-fractal (Fig.|8^) and the Apollonian network (Fig.|8j3). 
For both networks, one sees that coherence sets in for a > 1.5, 
and only beyond a certain threshold of the coupling strength. 

In particular, one observes the remarkable fact that coher- 
ence appears only in an intermediate coupling range, i.e. nei- 
ther to large nor to small values. This is in agreement with pre- 
vious works 1 18] concerning other systems of coupled chaotic 
oscillators, where one also observes that synchronized chaos 
requires that the coupling must be neither too weak nor too 
strong in order to avoid triggering spatial instabilities fz^. 
From Figs. Is); andlSjl one observes that all these transitions to 
coherence are of first-order. 
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For the pseudo-fractal, the upper threshold disappears when 
the heterogeneity a is further increased. However, for the 
Apollonian network the upper threshold not only persists but it 
shifts toward smaller and smaller coupling strengths when a is 
further increased. As far as we know, this is the first time that 
one observes such behavior, and it should be related to the ge- 
ometrical differences between both networks. In other words, 
since Apollonian is a very particular scale-free network, being 
the only one, studied so far which is embedded in Euclidean 
space, this particular feature seems to enable nontrivial syn- 
chronization behavior: stronger dominance in the coupling to 
the most connected nodes destroys coherence. 

IV. DISCUSSION AND CONCLUSIONS 

In this paper we studied fully synchronized solutions for 
three scale-free network topologies. The main conclusion is 
the following: in random scale-free networks synchronization 
of chaotic maps not only depends on the coupling strength but 
is mainly controlled by the outgoing connectivity k, which 
is a measure of cohesion in the networks. Because of that, 
one finds coherent solutions in random scale-free networks of 
fully chaotic logistic maps (a ~ 2) with outgoing connectiv- 
ity fc = 8 and homogeneous coupling, but not in determinis- 
tic scale-free networks, since they have rather small effective 
outgoing connectivity, namely fc = 2 for the pseudo-fractal 
network and fc = 3 for the Apollonian network. 

Therefore, although the exponent 7 of connection distri- 
butions in scale-free networks does not depend on the out- 
going connectivity we have shown that, in general, syn- 
chronization of chaotic maps in such coupling topologies is 
quite sensitive to it. Moreover, the transition to coherence is 
of first-order, indicating a similarity with other complex net- 
works In particular, the threshold values of the coupling 
strength obey a power-law, Eq. Q, as function of the outgo- 
ing connectivity. The exponent of this power-law depends on 
the nonlinearity a of the chaotic map, being almost constant 
below flc ~ 1.7 and decreasing linearly above it. Interest- 
ingly this value of Oc is in the vicinity of the bifurcation of 
the quadratic map where the period-3 window appears, and 
coincides with the appearance of other nontrivial behaviors in 
coupled map lattices with regular topologies, namely in the 
velocity distribution of traveling wave solutions L4ZI . 



The synchronization criterion was based here in the square 
mean standard deviation following previous studies llz^ . Of 
course, it could be possible to have numerically a zero stan- 
dard deviation ct^ ^ lO^'^" with a particular oscillator slightly 
nonsynchronized. However, in such case the deviation of the 
oscillator amplitude from the rest of the network would be of 
the order of 10"^^, a majorant of the precision for our crite- 
rion. Therefore, we believe that within this numerical preci- 
sion there are no spurious results. Further investigations could 
be done, implementing extensions of clustering criteria such 
as, e.g., that of Pikovsky et al. 1I43I1 . 

For deterministic scale-free networks with homogeneous 
coupling, the same value Uc indicates the threshold above 
which no coherent solutions are observed, independently of 
the coupling strength. Above Ac, coherence is observed only 
for heterogeneous coupling, namely for a > 1.5. However, 
for this range of values, we have also shown that coherence 
is also absent either for very small or for very large coupling 
strengths, due to spatial instabilities. Another particularly in- 
teresting result that still needs to be explained is that, for Apol- 
lonian networks, the coupling threshold beyond which coher- 
ence disappears gets smaller when the heterogeneity is further 
increased. This point is not observed for the pseudo-fractal 
network and may be due to the geometrical differences be- 
tween both deterministic networks. 

As a general property, we have shown that all transitions 
to coherence are of first-order. Furthermore, all results are 
robust not only against changes of the initial configurations 
of node amplitude but also, in random scale-free networks, 
against changes of the connection network. We also presented 
preliminary results indicating that in scale-free networks hubs 
play apparently no fundamental role in the dynamical collec- 
tive behavior, what remains to be further investigated. 
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